Raman cooling and heating of two trapped Ba + ions 
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We study cooling of the collective vibrational motion of two 138 Ba + ions confined in an electro- 
dynamic trap and irradiated by laser light close to the resonances S1/2-P1/2 (493 nm) and P1/2-D3/2 
(650 nm). The motional state of the ions is monitored by a spatially resolving photomultiplier. 
Depending on detuning and intensity of the cooling lasers, macroscopically different motional states 
corresponding to different ion temperatures are observed. We also derive the ions' temperature from 
detailed analytical calculations of laser cooling taking into account the Zeeman structure of the en- 
ergy levels involved. The observed motional states perfectly match the calculated temperatures. 
Significant heating is observed in the vicinity of the dark resonances of the Zeeman-split S1/2-D3/2 
Raman transitions. Here two-photon processes dominate the interaction between lasers and ions. 
Parameter regimes of laser light are identified that imply most efficient laser cooling. 
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I. INTRODUCTION 



In order to implement conditional quantum dynamics 
of qubits represented by the internal structure of indi- 
vidual ions in a Paul trap, a vibrational mode of the 
entire ion string is used. Cirac and Zoller Q have sug- 
gested a method for conditional dynamics that requires 
the ions to be cooled to their vibrational ground state. 
Here, the vibrational state represents a "bus" qubit that 
gets entangled with the ions' internal states by stimu- 
lated transitions, where the absorption (or emission) of 
a photon is coupled with the emission (or absorption) of 
a vibrational quantum. 

The recently proposed S0rensen-M0lmer scheme Q in- 
corporates two-photon transitions in order to directly 
entangle the internal states of the ions. Here, the vi- 
brational mode is only virtually excited. In this way 
the contribution of higher vibrational excitations can be 
made to destructively interfere. In this way, the scheme 
does not require cooling to the vibrational ground state, 
but admits population in higher vibrational levels. 

In general, the manipulation of the motion of trapped 
particles by light relies on the coupling of electronic and 
motional degrees of freedom via the scattering of pho- 
tons. Here, the important figure of merit is the Lamb- 
Dicke parameter (LDP) 77, its square given by the ratio 
of the recoil energy transfered when scattering a photon, 
and the energy of the vibrational quantum |1 . For small 
77, transition processes involving the creation or annihi- 
lation of n vibrational quanta (i.e. nth order sideband 
transitions) are suppressed by the factor rj 2n . Therefore, 
a small LDP allows only linear coupling (involving only 
single-phonon transitions) between electronic and vibra- 
tional degrees of freedom. This Jaynes-Cummings type of 
coupling |Q is prerequisite for many schemes suggested in 
quantum information processing including the mentioned 
implementations jl],|]]. 

Plcnio and Knight j^] have proposed to drive the in- 
teraction between electronic and vibrational degrees of 



freedom by tuning strong laser light to the resonance fre- 
quency, instead of the sideband. In this case, the light 
intensity imposes a specific light shift on the ions that se- 
lects the desired transition. By doing so, the time needed 
for conditional dynamics is reduced, as compared with 
the Cirac-Zoller scheme, thus reducing the effect of deco- 
herence. Here, too, the LDP plays a crucial role as the 
fidelity of such a fast quantum gate is / w 1 — \r} 2 . 

While the S0rensen-M0lmer scheme does not require 
ground-state cooling, it still demands low vibrational 
quantum numbers (i.e. r] 2 (n + 1) <C 1). Thus, laser cool- 
ing is essential for quantum information processing in ion 
traps, and it is desirable to have at hand robust cooling 
techniques that are suited for reaching low vibrational 
excitation at moderate technological expenditure on the 
side of the laser sources. In particular, cooling of the col- 
lective motion of several particles is prerequisite for im- 
plementing conditional quantum dynamics with trapped 
ions. 

It has been known for some time that the motion of 
trapped particles strongly depends on laser cooling. With 
a single trapped particle, distinct macroscopic motional 
states have been observed with different parameters of 
the cooling lasers [0 . With a cloud of ions, changes in the 
temperature induced by laser cooling cause phase transi- 
tions between a crystalline and an evaporated state [|],|| . 
These macroscopically distinct states indicate the specific 
temperature achieved by laser cooling. 

In this paper we investigate the macroscopic motional 
states of two trapped 138 Ba + under the action of laser 
cooling. In Sec. ||, the theory of laser cooling of the 
collective motion of two charged particles trapped in an 
anisotropic harmonic potential is outlined. If the par- 
ticles obey the Lamb-Dicke limit, the cooling rate and 
the achievable temperature are calculable without any 
free parameter. We fully take into account the multi- 
level structure of the involved ions, thus allowing for any 
type of cooling process, such as Doppler ||, sideband |lC| ] 
and Raman cooling |11|. In Sec. Ill we describe the ex- 
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perimental setup and present the results of the measure- 
ments. Phase transitions between macroscopically dis- 
tinct states are observed. The parameter values at which 
these transitions take place perfectly agree with the val- 
ues calculated from the theoretical model. In Sec. |^ 
we use the model to find optimized laser parameters for 
cooling the collective motion of the ions. 

II. THEORY 
A. Collective Motion 



the ions can move freely in a ring-shaped region in the 
(x, y) plane. However, at low temperature, the relative 
motion is restricted to the wells, i.e., the ions will appear 
at approximately fixed positions. The same argument 
can be applied to the z-direction: For even higher ener- 
gies E cx > AV Z the motional state changes from a ring to 
a sphere. All three motional states have been observed 
in our experiment. However, in the following we shall fo- 
cus on the transition between the localized and the ring 
state, as we are mainly interested in laser cooling in the 
low-temperature regime. 



We consider the collective motion of two ions of mass m 
in a three-dimensional (3D) harmonic trapping potential 



Krapfo) = f {u 2 x x] + uly) + U 2 Z Z 2 ) 



(1) 



which is the same for both ions (j = 1,2). Since each 
of the ions carries one positive elementary charge e, they 
repel each other by the Coulomb interaction 



Vc(ri,r 2 ) = 



1 



4ire ri - r 2 



(2) 



The potential is simplified by using center-of-mass and 
relative coordinates 



R 



ri + r 2 



and r = r 2 — ri . 



(3) 



With the definition of total and reduced masses M = 2m 
and jj, = m/2, one gets 



V(R) = f LjIx 2 + lo 2 Y 2 + u 2 z Z 2 



(4) 



for the center-of-mass motion, corresponding to the mo- 
tion of a single trapped particle with the original oscilla- 
tion frequencies, but twice the ions mass. 
The potential for the relative motion is 



V(r) = 1^ 2 xX 2 +u J y+u; 2 z z 2 ) + 



1 



47re r 



(5) 



If the trap frequencies slightly differ (as they do in our 
experiment, see Table this potential does not show 
spherical symmetry. With lo x < uj v < uj z holding, the 
potential has absolute minima on the x axis and saddle 
points on the y and z axis at 



r o, g =±\l 



47re \uJ\ q 



x,y,z 



(6) 



Figure |l] shows the potential in the (x. y) plane. It is 
clear that for high enough temperature, i.e. at a mean 
motional energy E ex comparable with the potential bar- 
rier at the saddle point on the y-axis 



AVy 



(J) 




FIG. 1. Effective double well potential V(x,y) for the 
relative motion of two trapped 138 Ba + ions (trap frequen- 
cies as given in Table |). The contour lines correspond to 
V = 2,4,6,8^eV. 



In the localized state the potential can be approxi- 
mated around the minimum at ro. x - Expanding up to 
second order in the displacement f = r — ro jX one gets 

V(r) = V(r , x ) + f {u 2 x x 2 + u 2 y y 2 + u 2 z ~z 2 ) (8) 

with local frequencies 

Co x = ViuJx , w 9 , 2 = Juj^ z ~ (9) 

(c.f. Table |) 



Mode 


X 


Y 


Z 


X 


y 


z 


i^/27rMHz 


1.0035 


1.0220 


1.0530 


1.7381 


0.1936 


0.3191 


f?493 


0.044 


0.044 


0.043 


0.034 


0.101 


0.079 


f?650 


0.034 


0.033 


0.033 


0.026 


0.077 


0.060 


f?1762 


0.012 


0.012 


0.012 


0.009 


0.028 


0.022 



TABLE I. Frequencies and Lamb-Dicke parameters for the 
different vibrational modes. Note that the geometrical factor 
(k;-e 9 ) is included in the LDPs (Eq. |l5|). In the present exper- 
iment the angle of incidence is 45° to all directions, yielding 
a factor of l/v3. 
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B. Laser Cooling of the collective motion 

The mechanical effect of irradiating ions with laser 
light leads to thermalization of the motional degrees of 
freedom. The corresponding cooling rate and steady- 
state temperature depend on the electronic level scheme 
and laser parameters; they are calculated in this section. 
The key to the theory presented here is the well-known 
Lamb-Dicke limit: As long as the ions are localized near 
the minima of the potential, their spatial extension is 
much smaller than the wavelength of the exciting laser 
light. In this case, the photon recoil operator e jki ' r , i.e. 
the coupling of the the internal dynamics to the ions mo- 
tion, can be expanded in the LDP. Since the LDP is small, 
dynamics due to this coupling is slow compared with the 
internal dynamics that may be adiabatically eliminated. 
This reduction leads to a master equation for the mo- 
tional degrees of freedom, where the involved transition 
rates depend on steady-state expectation values of inter- 
nal operators. This procedure and its application to a 
single Ba + ion has been described in [Q. 

We do not account for dipole-dipole interaction be- 
tween the two ions, since the separation of the two ions 
is large as compared to the optical wavelengths used 
(ro/A « 8.) Even though processes caused by this in- 
teraction possibly could be observed O], they are rare, 
compared with the single-ion scattering, and their influ- 
ence on the collective motion is negligible. Thus, we do 
not consider correlations between the two ions' internal 
dynamics that would lead to cooperative effects such as 
superradiance. We rather calculate the cooling effect that 
each of the two ions individually has on the common mo- 
tion. This effect is the same for both ions, since the foci 
of the lasers are wide compared with the distance of the 
ions (the waist of the Gaussian beams being w sa 10 ro), 
and therefore the laser intensities, and consequently laser 
cooling do not depend on the ions' positions. 

We shall now give a short recap of the theory presented 
in Ref. jl^] and apply its results to the collective motion 
of two ions in the current experiment. We start from the 
density operator g for the motional and internal degrees 
of freedom. The dynamics is described by the master 
equation 



Cg(t) 



(10) 



the index / runs over all atomic transitions involved. 2ji 
is the spontaneous emission rate, and af' are the cor- 
responding transition operators. Wj(£) is the angular 
distribution of spontaneous emission for this transition. 
The Hamiltonian is given by 

P 2 P 2 



H= 2M + V{R) + Yu + V{V) + £ H ^ + V »™ ■ 

^ 3=1,2 

(12) 



H nl _j is the Hamiltonian for the internal levels of the j-th 
ion. The ion's dipole interaction with the laser light can 
be written as 



V 



Dip ,j 



= > hfli e 



'1,3 



H.c. 



(13) 



where Sli is the Rabi frequency of the light driving the 
transition /. Summing over both ions, and rewriting the 
dipole interaction in terms of center-of-mass and relative 
coordinates [(Eq. g)], gives 

^ Dip = £ hCli e ik,R ( e-* ki r a+ + e* k ' r a+ ) + H.c. 
I 

(14) 

Note that the center-of-mass and relative motions are 
coupled via the interaction with the laser light. However, 
in the Lamb-Dicke limit this coupling can be neglected, 
as (k; • R) , and (k; • r) are small and the potential can be 
expanded. This expansion is characterized by the Lamb- 
Dicke parameters 



r)i, q = k; • e q yh/2muj q 



(15) 



for the center-of-mass modes and r\i q — k; • e q y/h/2muj q 
for the normal modes of the relative motion. In the ex- 
periment, all LDPs arc much smaller than unity, see Ta- 
ble Q. Now, we rewrite the coordinates in terms of the 
lowering and raising operators of the corresponding vi- 
brational eigenmodes 




where 



7 = 1,2 I ^ 



(ii) 



is the Liouvillean allowing fo r sp ontaneous emission. It 
has the usual Lindblad form H|, where dissipation acts 
on the internal (via a^) and on the motional (via e ±lk r ) 
degrees of freedom. The index j labels the two ions, and 



Since the expansion is around the potential minimum, 
there is a constant phase factor e ± ^ k, r °, that is incor- 
porated in the atomic transition operators a^-. Now, 
keeping only terms up to first order in (any) Lamb-Dicke 
parameter yields 



3=1,2 



V2 



(17) 
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where a x = a + + o~ and a v = (a + — cr~)/i have been 
used. The first term of Eq. ( |l7|) describes the driving of 
the internal dynamics by the interaction with the laser 
light. Since it is of zeroth order in any LDP it does not 
contain the mechanical effect of light, i.e. no laser cool- 
ing or heating. The second and third term of Eq. (|l7j ) 
(first order in a LDP) couple the internal dynamics to the 
center-of-mass and relative motions, respectively. The 
ion-dependent sign in the third line of Eq. ([l7|) arises 
from the opposite impact of the two ions' recoil on the 
relative motion. Since the vibrational modes decouple 
in the Lamb-Dickc limit, we can restrict ourselves to a 
ID description. Therefore, the index q characterizing the 
specific mode is omitted in what follows. 

Now, the fast internal dynamics of the ions is adia- 
batically eliminated and one obtains a master equation 
for the external degrees of freedom. This procedure fi- 
nally leads to the well known time evolution of the mean 
phonon number |p.5| 



j t (n) = -{A_-A + )(n)+A + 



with the cooling rate 



W = A- - A 4 



(18) 



(19) 



and (if the cooling rate is positive) the mean motional 
energy 



.Eex = hu>- 



A A 



A--A+ ■ 
The transition rates A± are given by 
A± = 2(ReS(=Fw) + D) 



(20) 



(21) 



The first constituent is the fluctuation spectrum of the 
electronic dipole coupling 

/>oo 

s(«) = rnifl PW / dte^Vf^KW), (22) 

where u) is the vibrational frequency of the considered 
mode. The expectation value is to be taken in the steady 
state of the internal dynamics that in turn is calculated 
from the master equation 



d , . 1 



H in + J2 m ^l > Qin(t) 



+ £o£>in(i) 



(23) 



derived from Eq. (|10|), by using adiabatic elimination, 
and keeping only the zeroth order in any LDP. Here, 

£oftn = 7i(2o-j~ Qi n af - o/crf g in - g- m a^a^) (24) 



describes the spontaneous emission in zeroth order. Via 
the driving by laser light, the steady state depends on 



the parameters (intensity, detuning, polarization) of the 
light field. 

The second contribution to the transition rates orig- 
inates from the spontaneously emitted photons as de- 
scribed by the Liouvillean L. So far, spontaneous emis- 
sion has only been considered up to the zeroth order in 
the LDP. The second order term is readily calculated and 
leads to the additional transition rate 



o 7 



(kj • e q 



(25) 



both from |n) to \n + 1) and vice versa. Here, the geo- 
metrical factor (k; • e q ) 2 in the LDP (eq. 15) is replaced 

by a; = J .. d£Wi(C)(, 2 that accounts for the angular dis- 
tribution of the light emission. 



III. EXPERIMENTAL OBSERVATION 

The two Ba + -ions are confined in a 1-mm-diameter 
Paul trap and irradiated by laser light at 493 nm for ex- 
citation of resonance scattering on the 2 Si/2- 2 Pi/2 tran- 
sition. This laser is detuned slightly below resonance for 
cooling the ions. The fluorescence signal is recorded by 
photon counting. A second laser at 650 nm prevents op- 
tical pumping into the 2 D 3 / 2 level. A static magnetic 
field defines the quantization direction and lifts the de- 
generacy of the magnetic sublevels. Its direction is at 
45 degrees with respect to the plane defined by the ring 
electrode. The directions of propagation and the polar- 
izations of both light beams are set perpendicular to the 
magnetic field. The power levels of the light fields are 
stabilized by electro-optic modulators. Precise detuning 
of the 650-nm laser is accomplished using a double-pass 
acousto-optic modulator. So far, the setup has been de- 
scribed elsewhere 16 1. In addition the fluorescent light 
is monitored by a spatially resolving photo multiplier re- 
vealing the position of the ions in the trap. 

The ions are dynamically trapped by superimposed 
electric dc and rf-fields |jl7) and represent particles inside 
a nearly spheroidal harmonic (quasi)potential. In order 
to determine the geometry of this (quasi)potential oscil- 
lations of a single trapped ion are excited by applying 
an additional rf field. Three sharp resonances are found 
at the frequencies given in Table [j], i.e. macroscopic mo- 
tion of the ion in the x, y and/or z direction of the trap 
is observed. The x-direction, within the radial plane, is 
found at approximately 45 degrees to the polarization of 
the lasers. While the nondegeneracy of uj z is obvious, the 
asymmetry of u x and oj v originates from a slight asym- 
metry of the ring electrode. The distance of the two ions 
in the potential minima, as calculated from Eq. (|J) using 
the actual values of \x and u) q is r — 3.70 /im and agrees 
very well with the experimental value 3.65 /im obtained 
from the spatially resolved measurement. 
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FIG. 2. Two trapped Ba ions show different motional 
states depending on laser parameters. Top: fluorescence of 
two trapped ions as a function of laser detuning, collected 
in 0.1s. Bottom: spatial distribution of the two ions at the 
detunings indicated above. 
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FIG. 3. Top: observed motional states for different detun- 
ings of the 650-nm light. The dots correspond to individual 
observations. Middle: mean motional energy in the y-mode 
calculated from theory. Bottom: cooling rate for the y-mode 
calculated from theory. 



Keeping all other parameters constant, the frequency 
of the 650-nm laser is scanned across the 2 D 3 / 2 - 2 Pi/2 
line. The intensity of the scattered light vs detuning of 
the 650-nm light, Agso = ^650 — wphD, is shown in the 
upper part of Fig. ^, where W650 is the laser frequency and 
^p^d is the atomic transition frequency. The structure 
of dark resonances is well known from the fluorescence 
of a single ion Jul and corresponds to four dark states, 
i.e. coherent superpositions of magnetic substates of 
the 2 Si/ 2 and 2 D 3 / 2 levels, respectively. From a fit of 
this excitation spectrum with eight-level optical Bloch 
equations, accounting for the Zeeman substructure, we 
obtain the Rabi frequencies of both lasers, the detuning 
A493 = W493 — ws„p of the 493 nm laser as well as the 
strength of the magnetic field. This is a standard proce- 
dure in our Ba + experiments. The parameters obtained 
from the spectrum shown in Fig. || arc Qs5o = 59.0 x 2n 
MHz, fZ 493 = 46.2 x 2n MHz, A 493 = -44.3 x 2tt MHz 
and B/h/iB = 7.7 x 2ir MHz, where /j,b is the Bohr mag- 
neton. 

Additional information is derived from monitoring the 
motion of the two ions. Depending on the detuning of 
the 650-nm laser, different motional states are observed. 
Whereas the ions appear at fixed positions if the red laser 
is tuned below the lowest or well above the highest dark 
resonance, they form a ring-shaped distribution in the 
(x, y) plane with certain intermediate detunings. 

The lower part of Fig. || shows the motional state at 
four different detunings, corresponding to the crosshairs 
in the spectrum. Note the clear distinction of the lo- 
calized state at Aq$q — —62.8 x 2ir MHz from the non- 
localized state only 2.6 x 2ir MHz away: Phase transitions 
between the two observed motional states appear at very 
definite laser frequencies. 



Obviously, the motional state of the ions depends on 
the ion temperature set by the laser cooling. In Fig. ||, the 
motional states observed at different detunings are shown 
together with the motional energies and cooling rates cal- 
culated from theory. The different motional states appear 
in distinct regions of detuning: The ions are observed in 
the fixed position for the detuning Agso below —61 x 2ir 
MHz. At this limiting value a sharp phase transition 
occurs, and the ions form a ring when —61 x 27r MHz 
< Agso < —56 x 2ir MHz. At the upper boundary, they 
are cooled again to the fixed-position state. Such a "hot" 
region appears with each dark resonance. At some val- 
ues of the detuning, e.g. at A 650 = — 19 x 2tt MHz, the 
motion fluctuated between the fixed and the ring-shaped 
state, thus remaining indefinite. 

The motional energy has been calculated from the 
model presented in the preceding section, using the ex- 
perimental values for the y-mode. This energy should 
be compared with the height of the potential barrier 
[(Eq. 0)], which for the given data is AV y = 7.15 fieV. 
Wherever the predicted temperature varies steeply upon 
detuning, a sharp transition is observed. This is the case, 
when the corresponding cooling rate changes from posi- 
tive to negative, i.e. when the lasers actually heat the ion. 
Laser heating appears at negative detunings of both lasers 
when Raman transitions dominate the interaction. In the 
vicinity of the dark resonances, the internal state con- 
sists of a coherent superposition of the S1/2 and the D 3 / 2 
states, forming an effective two level system. Now, with 
both detunings being negative, and Agso > A 4 g 3 + A zee , 
the effective detuning of this Raman transition is posi- 
tive, and causes heating of the motion. Here, A zec is the 
difference in Zeeman shifts of the involved sublevcls that 
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defines the resonance condition of the considered dark 
state. 

Note that the transition to the ring state, at positive 
detunings from the Raman resonances, appears where the 
expected energy is still considerable lower than the poten- 
tial barrier. However, this does not contradict the model, 
which assumes the ions close to the potential minima. As 
soon as they are smeared out over a region of considerable 
size, the harmonic approximation around the minimum 
no longer holds. The cooling gets less effective for the 
broadened velocity distribution. This happens when the 
thermal energy is substantially lower than AV y . Still, 
the observed phase transitions agree very well with the 
results of the model. 

At low speed of scanning, for example, dA^/dt ss 
1 x 2ir MHz/s, we observed hysteresis effects: The mo- 
tional state upon up-tuning the laser frequency differed 
from the state observed upon down-tuning. However, 
this effect vanished when the scanning rate was de- 
creased. Scans at a very low speed (dA e50 /dt < 0.1 x 2n 
MHz/s) showed that the motional state solely depends 
on the laser parameters. Therefore, a steady state exists, 
which is approached in a finite motional relaxation time. 

The final temperature is, in principle, not only set by 
the laser cooling, but it is subject to additional heat- 
ing mechanisms, such as fluctuations in the trapping 
electric field. However, the calculated transition rates 



A, 



10 5 



10 6 s 1 equal the number of scattered 



phonons per unit time, and thus exceed the rate of par- 
asitic heating (ps 10 . . . 10 4 s _1 ) observed in ion trap ex- 



periments |19|. Therefore, this additional heating can be 
neglected as long as the ions are kept laser-cooled. 



IV. OPTIMIZATION OF LASER COOLING 

The discussed model is now employed to numerically 
investigate the influence of various parameters on the 
laser cooling, and to find optimized conditions for ef- 
fective and robust cooling. We concentrate on the x 
(stretch)mode, since it is less affected by perturbations 
from stray fields of the trap |2(J and therefore is suitable 
for serving as a bus mode between the ions. A promis- 
ing candidate for implementing qubits in trapped Ba + 
is the Si/ 2 -D 5/ 2 (1762 nm) transition The LDP 

corresponding to this infrared transition is smaller than 
those of other transitions (see Table |) , and the S0rcnsen- 
M0lmer criterion r) 2 ((n) + 1) <C 1 yields the condition 
(ri) <C 10 4 that is easily fulfilled in practice. In fact, a 
mean phonon number (n) < 1 can be reached: This is 
shown in the middle part of Fig. |], where (n) is plotted vs 
the detuning of the 650- nm light. The lowest vibrational 
excitation was found by randomly searching the 4D pa- 
rameter space for the lowest achievable (n). Then, three 
laser parameters were kept constant, while the red-laser 
detuning was varied, in order to generate the graph. The 
minimal temperature is achieved around Agso ~ 6 x 2ir 



MHz. Note that the lowest value (n) — 0.9 indicated 
by the vertical dotted line does not match a particularly 
large cooling rate, as is shown in the lower part of Fig. |^. 
Looking at the expected fluorescence signal (upper part 
of the Figure), which corresponds to the population of 
the P1/2 level, we note that the optimum detuning is on 
the red side of the dark resonance at A^q = 4.9 x 2ir 
MHz, as is expected for Raman cooling. The dark res- 
onance is power broadened by the high intensity at 650 
nm. On the other hand, we find low temperatures only 
for low intensity of the 493-nm light. This fact indicates 
that scattering of the blue light is the dominant heating 
mechanism. We do not claim to have found the absolute 
minimum of achievable temperature: Other sets of laser 
parameters may yield similarly low temperatures, and 
also the applied magnetic field to some degree influences 
the achievable temperature. 




A 650 (2^MHz°) 

FIG. 4. Fluorescence, mean number of vibrational quanta 
in the i-mode, and cooling rate vs laser detuning. The 
vertical dotted line indicates the minimal (n), laser param- 
eters: Sl 6 50 = 89.6 x 2tt MHz, f2 4 g 3 = 2.1 x 2vr MHz, 
A493 = —0.21 x 2-7T MHz. Other parameters as in Sec. HI 



We have restricted the parameter search to the range of 
negative detuning. Even lower energies may be reached 
with both lasers detuned to the blue side of the atomic 
resonance and the parameters adjusted to yield so-called 
electromagnetically induced transparency cooling [B2[ . 
Here, the dynamic Stark contribution to the nonlinear 
susceptibility dominates the interaction of ions and light 
p3| . However, Doppler precooling is still necessary and 
the lasers must either be switched from negative to pos- 
itive detuning, or else additional light sources are re- 
quired. Moreover, the intensities of the lasers must be 
precisely controlled, since for this type of cooling to work, 
the dynamic Stark shift should match the vibrational fre- 
quency particularly well. However, cooling mechanisms 
are desirable, that are more robust against variations of 
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intensity and detuning of the laser light, and thus are 
easier to implement experimentally. 

The robustness of the cooling is affected by the sharp 
boundaries between parameter regimes of very low and 
high temperature, because in these regions any drift of 
the lasers' frequencies may cause drastic deterioration of 
the cooling. In order to find "safe" parameters we re- 
peated the search, incorporating a tentative drift in the 
following way: For a given point in the (A493, Ag5o) pa- 
rameter plane, (n) was calculated at this point and at 
eight additional points, where one or both of the detun- 
ings differ by ±1 x 2tt MHz, and then averaged over the 
(n) values at all nine points. The best robust parameters 
found differ from the ones used to calculate Fig. || and 
the achievable ion temperature is not as low as in the 
ideal situation, but one still finds (n) = 3.3. The result 
is shown in Fig. |^, where the mean phonon number is 
shown as a function of A650 (solid line) . The dotted line 
shows, for comparison, the vibrational excitation that is 
achieved with the same parameters, when the averaging 
over nine points is not carried out. The regions of safe 
cooling are narrower, since detuning close to a "bound- 
ary" between cooling and heating is penalized. It seems 
not advisable to use a rather narrow region of cooling (as, 
for example, close to A650 ~ — 7 x 2ir MHz), because spu- 
rious laser drift might easily impose unexpected heating. 
However, some regions [e.g. A^o w (—22 to —17) x 2tt 
MHz and A650 ~ (—5 to 5) x 2tt MHz] are broad enough 
to allow for some laser drift with satisfactory cooling. 
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FIG. 5. Mean number of vibrational quanta for fixed 
laser parameters (dotted) and averaged over nine combina- 
tions of both detunings (see text) (solid). Laser param- 
eters: Sl 6 50 = 59.5 x 2tt MHz, 0,493 = 3.8 x 2tt MHz, 



A493 = —10.9 x 2-7T MHz. Other parameters as in Sec. [II 



V. CONCLUSION 

We have observed the effect of laser cooling and heat- 
ing on the collective motion of two trapped Ba + ions 
by monitoring the macroscopic motional states of the 
ions. These states indicate the vibrational excitation, 
or temperature of the ions. In this way, Raman heating 
is demonstrated that dominates on the blue side of the 
dark resonances (i.e. Agso > A 4 g3 + Az co ). 

The observed phase transitions are a remarkable 
macroscopic manifestation of dark states, i.e. of the co- 
herent quantum superposition of atomic levels. Raman 



cooling on the red side of the dark resonance is implied, 
though it is not directly observed with this method, since 
no macroscopic change of the motional state shows up. 

The good agreement of the calculated excitation with 
the experimental observations shows that the outlined 
theory describes well the laser cooling of the collective 
motion of the ions. The predicted cooling rates much ex- 
ceed spurious heating rates of the vibrational motion, for 
example, from fluctuations of the trapping electric field. 
Thus, such heating is negligible as long as the cooling 
light is on. 

The model comes with no free parameters; however, 
it is applicable only if the ions are localized near their 
potential wells. Nevertheless, this is not a strong restric- 
tion since this regime is supposed to be used in any ex- 
periment implementing conditional quantum dynamics. 
This allows us to identify parameter regimes that im- 
ply effective and robust laser cooling. We find that the 
collective motion can become cooled close to the ground 
state, if the dark resonances are not much broadened by 
laser power or fluctuations. But even when fluctuations 
in laser parameters are allowed for, low levels of temper- 
ature emerge that are suitable for quantum information 
processing using the scheme for conditional quantum dy- 
namics suggested by S0rensen and M0lmer M. We note 
that the level scheme of Ba + has the same structure as, 
for example, Ca + p3j, a type of ion often suggested, and 
the results obtained here are applicable, with appropriate 
modifications, to the cooling of these ions, too. 
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